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Abstract. In this paper, we study the motion of the free surface of a body 
of fluid over a variable bottom, in a long wave asymptotic regime. We assume 
that the bottom of the fluid region can be described by a stationary random 
process I3{x, u>) whose variations take place on short length scales and which 
are decorrelated on the length scale of the long waves. This is a question 
of homogenization theory in the scaling regime for the Boussinesq and KdV 
equations. 

The analysis is performed from the point of view of perturbation theory 
for Hamiltonian PDEs with a small parameter, in the context of which we 
perform a careful analysis of the distributional convergence of stationary mix- 
ing random processes. We show in particular that the problem does not fully 
homogenize, and that the random effects are as important as dispersive and 
nonlinear phenomena in the scaling regime that is studied. Our principal re- 
sult is the derivation of effective equations for surface water waves in the long 
wave small amplitude regime, and a consistency analysis of these equations, 
which are not necessarily Hamiltonian PDEs. In this analysis we compute the 
effects of random modulation of solutions, and give an explicit expression for 
the scattered component of the solution due to waves interacting with the ran- 
dom bottom. We show that the resulting influence of the random topography 
is expressed in terms of a canonical process, which is equivalent to a white 
noise through Donsker's invariance principle, with one free parameter being 
the variance of the random process f3. This work is a reappraisal of the paper 
by Rosales & Papanicolaou |24l and its extension to general stationary mixing 
processes. 



1. Introduction 

The problem of surface water waves over an uneven bottom is a classical prob- 
lem of fluid mechanics, and it is relevant to coastal engineering and ocean wave 
dynamics. In this paper, we investigate how the presence of bottom topography 
affects the equations describing the limit of solutions in the long wave regime. We 
assume that the bottom is modeled by a stationary random process which is mix- 
ing, whose variations and whose correlation length manifest themselves on length 
scales that are short compared to the scale of the surface waves. In a previous 
work [9], we addressed the long wave limit of surface waves over a bottom which 
has periodic variations over short scales, in which we proved that the problem fully 
homogenizes. That is to say, the free surface motion can be described by a par- 
tial differential equation with constant effective coefficients, where the dependency 
over short scales is manifested by coefficients which are ensemble averages. Here 
in contrast, we show that random, realization- dependent effects are retained in the 
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description of the solution. The latter paper and the present one are reappraisals 
and extensions of an earlier work by Rosales & Papanicolaou [24] who address the 
problem through different methods. 

Our approach uses a formulation in terms of perturbation theory for Hamiltonian 
partial differential equations, coupled with a detailed analysis of stationary ergodic 
processes which have mixing properties and which are considered as tempered dis- 
tributions. As a first result we give an appropriate form of the Boussinesq equations. 
Secondly, following a series of changes of variables, we derive a system of coupled 
KdV-like equations for the two components of the solution; these describe a wave 
propagating predominantly to the right, and a 'small' scattered wave propagating 
to the left. We then extract a limiting system of two effective equations through a 
consistency analysis. Specifically, we solve the effective system, which is composed 
of an equation similar to the KdV for the wave propagating to the right with a 
random component to its velocity, and a scattered wave propagating to the left. 
We give explicit formulas for the dominant contributions and the first corrections 
to this solution, quantifying the effects of the random modulation of position and 
amplitude. From these expressions, we compute a posteriori all the terms that have 
been neglected in the effective system, and prove that they are indeed of higher or- 
der. This evaluation relies on scale separation lemmas, which in turn follow from 
Donsker's invariance principle. Our analysis improves upon [24j in several ways. 
In particular we identify the canonical limiting distributions which contribute to 
the random asymptotic behavior of solutions, we quantify both random phase and 
random amplitude variations of solutions, and in addition, we extend the long- wave 
analysis over random topography to general stationary mixing processes. 

The asymptotic system of equations that results from this analysis consists of 
a KdV equation with an additional linear term, and a transport equation for the 
scattered component driven by an inhomogeneous forcing term. The additional 
nonzero linear term, which either stabilizes or destabilizes solutions depending upon 
the sign of its coefficient, in turn depends on the statistics of the bottom variations. 
The presence of this term is the consequence of a subtle calculation, and to our 
knowledge, it has not been previously observed. In case these statistics are spatially 
reversible, the relevant coefficient vanishes and the equation reduces to the usual 
KdV. 

There has been a lot of interest in wave motion in basins with non constant 
bathymetry, due to its hydrodynamic importance. Recent references to the theory 
of linear waves include the papers of Nachbin (1995) fT^, S0lna & Papanicolaou 
(2000) 26 , Nachbin & S0lna (2003) ,21J which discuss the theory of linear trans- 
port in a random medium. The earlier work of Howe (1971) |15j and the paper of 
Rosales & Papanicolaou (1983) [24] give an asymptotic analysis of nonlinear equa- 
tions of water waves. Nonlinear problems over variable topography are addressed 
in Nachbin (2003) [ID] and Artiles & Nachbin (2004) [2]. More recent contributions 
which take into account the combined effect of randomness and nonlinearity include 
the series of papers by Mei & Hancock (2003) JJj and Grataloup and Mei (2003) 
[l4] on the modulational scaling regime, and its extensions to the three dimen- 
sional case in Pihl, Mei & Hancock (2002) [23] . This work focuses on the temporal 
behavior of ensemble averages of solutions, giving the result that they satisfy a 
nonlinear Schrodinger equation with an additional dissipative term. The analog of 
this picture in the long wave scaling regime appears in Mei & Li (2004) [18] , where 
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the bottom is assumed random but varies on the same spatial scale as the surface 
waves. 

There is a history of rigorous analysis of the initial value problem and limiting 
equations in the long wave asymptotic regime of the water wave problem. Most 
of this work concerns the case of fluid domains with a flat bottom. The papers 
that address the KdV limit include Kano & Nishida (1986) [16^, Craig (1985) [7], 
Schneider & Wayne (2000) [25,, Wright (2005) ^ and Bona, Colin & Lannes 
(2005) 0. A recent paper which addresses specifically the Boussinesq scaling limit 
of the problem on a rigorous basis, and categorizes the well-posed possible limits 
is Bona, Chen & Saut (2002) [4]. There has been several papers giving a rigorous 
analysis of the initial value problem of water waves over a variable bottom, including 
Yosihara (1983) f29| on the two-dimensional problem and Alvarez-Samaniego & 
Lannes (2006) [1^ on the two and three-dimensional problems, and a recent paper 
by Chazel (2007) [B]. The paper [T] considers the issue of convergence in various 
scaling regimes governed by long wave models. These results are in the context 
of a deterministic problem, with a small amplitude bottom perturbation, varying 
spatially on the same scale as the waves in the surface. As far as we know, there are 
no current rigorous analytic results for the KdV or Boussinesq scaling regimes in 
which the bottom variations occur on a short length scale, and are averaged under 
the nonlinear evolution of water waves. 

The paper is organized as follows. Section 2 describes the problem of water 
waves in its Hamiltonian form, the Dirichlet-Neumann operator in the presence of 
a variable bottom, and the spatial scaling regime appropriate for the long wave 
problem. Section 3 presents the setting of stationary ergodic and mixing processes 
in which we work, and gives the relevant scale separation lemmas. This is the 
key of the paper. It furthermore gives an analysis of the natural regularization of 
characteristic coordinates that are applied to the KdV scaling limit. The Boussinesq 
regime is presented in Section 4, while the more detailed KdV regime is taken upon 
in Section 5. The main issue of this analysis is that the scattering of waves by the 
bottom variations is strong and it must be shown that the standard KdV Ansatz of 
unidirectional propagation remains valid despite this. The consistency analysis of 
the resulting asymptotic system of equation is the most detailed part of this paper. 
Finally, Section 6 presents some remarks on the process of ensemble averaging. 



2.1. Hamilton equations. The time-dependent fluid domain consists of the re- 
gion S{l3,ri) = {{x,y) e M""^ xR : -h + P{x) < y < vix,t)}, in which the fluid 
velocity is represented by the gradient of a velocity potential. 



The dependent variable 77(2;, t) denotes the surface elevation, and f3{x) denotes 
the variation of the bottom of the fluid domain from its mean value. The bottom 
variations are chosen from a statistical ensemble P), which is indicated by the 

notation (3 = (3{x,uj). The details of the ensemble and the associated probabilistic 
properties are described in Section 13.11 

On the bottom boundary {y ~ —h+l3{x)}, the velocity potential obeys Neumann 
boundary conditions 



2. Hamiltonian formulation 



u = Vif , 



Aip = 0. 



(2.1) 



V^-iV(/3) =0 



(2.2) 
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where N{P) = (1 + \dxP\'^)^^^^idxP, —1) is the exterior unit normal. 

The top boundary conditions are the usual kinematic and Bernoulli conditions 
imposed on {(a;, y) : y — ri{x, t)}, namely 

dti] = dyf ~ dxf] ■ d^if, dtif = -grj- ^\\/ipf . (2.3) 

The asymptotic analysis in this paper is initiated from the point of view of the 
perturbation theory of a Hamiltonian system with respect to a small parameter. 
For this purpose we describe the water wave problem as a Hamiltonian system 
with infinitely many degrees of freedom. In [30j . Zakharov poses the equations 
of evolution (12. ip (|2.2p (|2.3p in the form of a Hamiltonian system in the canonical 
variables {r]{x),^{x)) where one defines £,{x) = ip{x,ri{x)), the boundary values of 
the velocity potential on the free surface. The evolution equations take the classical 
form 

a(;)^(° p.^) 



with the Hamiltonian functional given by the expression of the total energy 



H ^ I I -\Vip{x,y)\'dydx+ n^rf{x)dx 

-h+f3{x) ^ J ^ 

ie(x)G(/3, 7j)ax) dx + l^rj^ {x) dx . (2.5) 

The Dirichlet-Neumann operator G{f3, if) is the singular integral operator with 
which one expresses the normal derivative of the velocity potential on the free 
surface. It is a function of the boundary values S,{x) and of the domain itself, as 
parameterized by /3(a;) and 'rj{x) , which define respectively the lower and the upper 
boundaries of the fluid domain S{l3^rj). That is, let ip{x,y) satisfy the boundary 
value problem 

A(p = in S'(/3,7/) , (2.6) 
\/ip-N{(3) = on the bottom boundary {y^—h + P{x)}, 
ip{x,r]{x)) = ^(x) on the free surface {y = 77(3;)} . 

The Dirichlet-Neumann operator is expressed as follows 

G(/3, 7y)e(x) = \/^{x, 7j{x)) ■ N{rj){l + \d,rfy/^, (2.7) 

where N{r]) is the exterior unit normal on the free surface. It is clearly a linear 
operator in ^ and it is self-adjoint with this normalization. However it is nonlinear 
with explicitly nonlocal behavior in f3{x) and i^^x). The form of this operator, and 
its description in terms of /3 and t] are given in the next section. 

2.2. Description of G{/3,ri). We now restrict consideration to the dimension n — 
2. In the undisturbed case in which the bottom is flat, the solution is formally 
given by a Fourier multiplier operator in the j:-variable. Using the notation that 
dx = iD; 

'^^ J J cosh(fc/i) ' cosh.{hD) ' ^ ' 
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When the bottom topography is nontrivial, as represented by {y — ~h + /3(x)}, 
the expression ()2.8p is modified by adding a second term in order that the solution 
satisfies the bottom boundary conditions 

^ '°tish(t^)^^ ^^^^ + sinHyD){Lm){x) . (2.9) 

The first term in (j2.9p satisfies the homogeneous Neumann condition at y = — ft, 
while the second term satisfies the homogeneous Dirichlet condition at ?/ = 0. The 
operator L{(3) in the second term acts on the boundary data £,{x) given on the free 
surface. In |9] we analyzed i(/3) in a nonperturbative case, where \(3\c^ ~ C'(l)- 
Here we are restricted to a perturbative regime, where we describe the expansion 
of the operator G(/3, 77) for small but arbitrary perturbations rj{x) of the surface, 
and small bottom variations l3{x). 

At order and 0{t]), one gets G(°) ^ D tanh{hD)+DL{/3) and G^^^ = Dr/D- 
G(°)r?G(°). At higher order, one finds the same recursion formula for G^'-* as for the 
case of a fiat bottom [10] except that the role of the operator Gq = D tanh(/iD) is 
now replaced by G^°^ . 

Since we allow bottom perturbations to be of order 0(e), we will use a recursion 
formula given in [9] for L(/3) in powers of /?. 

L(/3) = Li(/3) + L2(/3) + .... (2.10) 

with the first terms being 

Li(/3) = -sech(fti:>)/3sech(/ii:>)i:> (2.11) 

L2(/3) = sech{hD)PDsmh{hD)Li 

= -sech{hD)(3Dtaiih{hD)f3Dsech{hD). (2.12) 

General formulas are presented in [9j together with a Taylor expansion of the 
Dirichlet-Neumann operator G(/3, rj) in powers of both f3 and 77. In the analysis 
of the present paper, we will need only the terms up to second order in /3. 
The Hamiltonian is thus expanded in powers of 77 and /3 in the form 

Hiv,C,f3) {^Dis.iMhD)^ + gif)dx 

-iy" £_Dsech{hD)(3Dsech{hD)^dx 

+ ^ J ^{DtjD - Dta.nh{hD)T]Dtanh{hD))S,dx 

~\ J ^i^^'^^HhD)f3Dtanh{hD)l3Dsech{hD))£_dx 

+o{p''e) + oivpe) + o{rfe) ■ (2.13) 

By integration by parts, 

H{r].X\(i) = - / {£,D timhQiD)^ + gff)dx - - [ f3\Dsech{hD)^\'^dx 



-i J ^{DrjD - D tanh{hD)i]D tanh{hD))^dx 



J {Daech{hD)()PDtanh{hD)PDsech{hD)^dx 

+0{/3'e) + OivPf) + 0{ife) , (2.14) 
which is the starting point for our asymptotic expansion. 
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2.3. Spatial scaling and the scaled Hamiltonian. We consider the case in 
wliich tlie bottom varies on a short length scale, that is (3 — P{x,oj) is a random 
process, of zero mean value that satisfies ergodicity and mixing properties which 
will be detailed below. 

The fundamental long wave scaling for the problem of surface water waves re- 
tains a balance between linear dispersive and nonlinear effects in the dynamics 
of the surface evolution. The scaling that anticipates this balance is through the 
transformation 

X = ex, £,{x)=el{X), r]{x) = e^fi{X). (2.15) 

As for the bottom, we assume its variations are of order 0{e), which are much 
larger that the variations of the surface elevation, namely 

I3{x,uj) = £j3{x,uj). (2.16) 

We assume that (3 is bounded in for almost every realization lo € Vl. 

In order to get the scaled Hamiltonian, we need to examine the asymptotic 
expansion of the Dirichlet- Neumann operator G(/3, 77) in a multiple scale regime. We 
recall how formally a pseudo-differential operator acts on a multiple scale function 
f{x,X) where X — ex (see [H] for details). In particular let m(D) be a Fourier 
multiplier operator acting on a function /, defined as 

{m{D)f){x) = i- y e''^^--y^m{k)f{y)dydk. (2.17) 

When m{D) acts on a multiple scale function f{x,X) with X = ex, D is replaced 
by Dx + eDx and 

m{D)f{x,X) = /e^M-^.)(^!!^Lie.^^)/(y,x)dydfc 

= m{D,)f + em'{D,)Dxf + ■■■ (2.18) 
Applying this to the scaled Hamiltonian, we get 

H{fi, h A f) = Y y {KDU + gfi^)dX (2.19) 

"Y j m\Dx^ech{ehDx)~^\^dX J ^{D x^lD xl - ^ DU)dX 

-y / {Dxscc\i{ehDx)0 
(3{x){D^ + eDx) t&nh{h{D^ + eDx))i3{x)Dxsech.{ehDx)i dX 

For simplicity of notation, we now drop the tildes over /?, 77, ^. Expanding the 
operator sech{ehDx) in the second term in (|2.19p gives 



J I3{x) \DxsecHehDx)^\^dX = J f3{j) 



Dxil-le^h'Dj,)C 



dX . (2.20) 



The last term of (|2.19p is a little more complicated but is calculated in the same 
manner. Expanding {Dx + eDx) ta.nh{h{Dx + eDx)) we get 

{Dx + eDx) tanh{h{Dx + eDx)) = D^ tanh{hDx) + 0{e) . (2.21) 
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Finally, 

/ Dxsech{ehDx)S. [f3{x){D^ + eDx) tanli(/i(L>^ + eDx))P{x)Dxsec\i{£hDx)£] dX 



= J Dxsech{ehDx)^[l3{x)D.^ tSinh{hD.^)f3{x)]Dxsech{ehDx)CdX + 0{e) 
= J[(3ix)D,tsiiih{hD,)Pix)] \Dx^\^dX + 0{s) . 

(2.22) 

Putting all these terms together: 



Hi-n, ^;/3,e)^^ J [(/i - eP{x) - e^P{x)D^ t'dnh{hD^)/3{x)^ \Dx^\^ + gif 
+ ^{^DxvDx^ - ^^D%d dX + o{e^) . 



(2.23) 



3. HOMOGENIZATION AND SCALE SEPARATION 



The purpose of this section is to understand the asymptotic behavior of integrals 
of the form 

f+oo -Y 

7(-)/(X)dX:=Z,(7,/) , (3.1) 

where f{X) comes from expressions which involve the physical variables which 
depend only upon large spatial scales, and where 7(3;) = 7(0;; w) is a stationary 
ergodic process taken from the statistical ensemble VL from which our realizations 
of the bottom are sampled. Principle examples of such integral expressions in the 
Hamiltonian for water waves are 

/+00 v" 
P{--Lo)\Dx£.{X)?dX (3.2) 
-00 £ 

as well as 

{l3D,i&nh{hD,)l3){-)\DxaX)? dX . (3.3) 

In our previous work [9] , expressions of this form are analyzed under the hypothesis 
that 13 was a periodic function of x. In the present paper, we are concerned with 
the case in which the bottom variations I3{x^ lo) are decorrelated over large spatial 
scales, which is quantified with a mixing condition on fJ. 

3.1. Stationary ergodic processes and mixing. We take our statistical ensem- 
ble of random bottom variations of the fluid domain to be modeled by a stationary 
ergodic process which will possess some properties of mixing. Mathematically, given 
a probability space (f2, M, P) equipped with a group of P-measure preserving trans- 
lations {tj, : y S M}, and a function G : ^ R, then a stationary process 7 is given 
by 7(0;; cj) :— G{txUj). The notation for the probability of a set A G A4 is P(^), 
and integrals of functions F over this probability space are denoted by 

FdP^E{F). (3.4) 

o 

We further require that the measure be ergodic with respect to {rj,}j,gii, meaning 
that for any bounded measurable function F : R, then for P-almost every 
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realization a;, 

lim J I F{TyLo)dy^E{F) . (3.5) 

For our purposes, we would like to take := C(E) the space of bounded con- 
tinuous functions, for which the one-parameter group of translations is just that, 
('^y7)(') — li' + y)j fo'' 2/ S However it turns out that our sample space C(K) 
must be enlarged to a subset of the space of tempered distributions S' , as the pro- 
cess of taking limits invokes Donsker's invariance principle, and the support of our 
limiting measures is on distributions corresponding to one (or several) derivatives 
of Brownian motion. The modeling of a random bottom will require properties 
of asymptotic independence of typical realizations with respect to the probability 
measure {Ai, P), specifically that the translations {Ty}y^]R exhibit a mixing prop- 
erty with respect to it. There are several notions of mixing in the literature [13j . 
For simplicity, we adopt the notion of uniform strong mixing (called a-mixing), 
although weaker conditions would also work in our setting. The stationary pro- 
cess defines a natural filtration on the probability space given by the cr-algebras 
tW" = (7(7(2/, w) ■ l£ y l£ u). The notion of a-mixing is that there is a bounded 
function a{y) for which a{y) ^ as y — > 00 such that for any two sets A E Ai^ 
and B e M°_^ then 

\PiAnTyiB))~PiA)PiB)\<aiy) . (3.6) 

Note that mixing implies the process is ergodic. So that Donsker's invariance prin- 
ciple will extend to this mixing process [22], we require that a{y) = 0{l/y\og{y)) 
for y I—!- +00 as well as 

POO 

/ a{y) dy < +00 . (3.7) 
Jo 

The integral (|3.3p involves a nonlocal expression in the bottom variations (3{x), 
implying that the random processes we are led to analyse will never be perfectly 
decorrelated under any finite translation. Indeed, the spatial decay of the kernel of 
the operator Dtaiih(hD) implies a lower bound on a{y) of the form 

aiy) > e-^^y , 

even for statistics of the actual realizations of the bottom variations [3{x,uj) which 
are fully decorrelated under sufficiently large finite translations \y\ > R. 
For the zero mean process 7, define the covariance function pj to be 

p^iy) E(7(0; Lu)jiy; lu)) = £(7(6; c^)t^7(0; u;)) , (3.8) 

which is an even function of y ([12] page 123, or [3], page 178). The variance fi^ is 
given by the expression 

/•oo 

cr^ := 2 / p-t{y)dy . 
Jo 

The integral exists because of the hypothesis of mixing of the underlying process. 
The variance can take on any value in [0,-|-oo), and we are principally concerned 
with the situation in which cr^ > 0. To this end we note the following fact. 

Lemma 3.1. When the process (3{x,uj) = dxj{x,uj), for 7(2:) G , a zero-mean, 
stationary process with the above mixing properties, then 
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Proof. By definition, 

r+oo r+oo 

al = 2 E{mPiy))dy = 2 EWx)P{x + y))dy (3.9) 

Jo Jo 

r+oo /• + 00 

= 2/ E{dj:^{x)dxjix + y))dy ^ 2 E{d^j{x)dyj{x + y))dy 

Jo Jo 

r+oo 

= 2 dyE{d,;j{x)j{x + y))dy . 

Jo 

Therefore by integrating, 

aj = -2E(a,7(x)7(a;)) + 2E(a,7(x)7(a; + y))\y=+oc = --E{d,j\x)) , 
because tlie process is mixing. Using tlie hypothesis of ergodicity, 

E(a.7'(^)) = lim ^ f d.,^^{x)dx = hm b-i\x)) ^ = . (3.10) 

T^co 1 Jq T^oo 1 x=0 

Thus the most interesting processes are those which are not derived from deriva- 
tives of another stationary process; this fact will be reflected in our analysis of the 
asymptotics of the integrals (|3.ip in the next section. □ 



3.2. Scale separation. The asymptotic analysis of Hamiltonians or partial differ- 
ential equations which involve random coefficients needs to establish a clear criterion 
with which to characterize terms by their order parameter. In our present analysis, 
we view each term as a tempered distribution in space and time, namely in iS'(M^). 
We consider a term a(X, t; e) to be of order 0{e^') if for any Schwartz class test func- 
tion (p{X, t) the limit lim^^o £ / o^i^, t', £)v(^i d^dt exists. In this context, the 
terms of a partial differential equation with random coefficients represent random 
ensembles of tempered distributions, say {a(X,t;uj,e) : oj E fl} C 5'(]R^), which we 
state to be of order 0{£'^) if for any test function (p{X,t) G 5(R^) the probability 
measures dP^ of / a{X,t; e,uj)(p{X,t) dXdt converges weakly to some dPo. In 
this section we discuss the behavior of such terms in the form 

/ 7(f , t; Lj)v{X, t)(p{X, t) dXdt 

J 7i(f , t; c^)72(^, t; cu)v{X, t)^{X, t) dXdt ^^'^^^ 

where 7 is a stationary mixing process, w is a solution to one of the several differ- 
ential equations under discussion, and Lp plays the role of a test function. 

Lemma 3.2. For ^{x\ljj) a stationary ergodic process and for f{X) G L"'^(M), then 
for P-a.e. realization uj, 

f{X)j{—,uj) dX = E(7) / f{X) dX + 0(1) . (3.12) 

-co £ J — 00 

Proof. For a Schwartz class function / we have 

+ 00 y r + co I 

f{X)^{--Lo)dX^e f{X) — { i{s-u:)ds)dX 



£ J-00 dX Jq 



+00 ^ 



Xf'(X)— I 'y{s;uj)dsdX. (3.13) 
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As £ — > 0, combining Birkhoff ergodic theorem 

^^\(s;c.)ds^E(7) (3.14) 
with the dominated convergence theorem leads to 

r+oc Y r+co 

f{X)j{—,Lu)dX ^-Ei^) Xf{X)dX (3.15) 

and finaUy (|3.12[) . In fact it suffices that / e L^(R) for the result to hold. □ 



The immediate application of the lemma is to the integrals (|3.2p ()3.3p , at least to 
the order implied by Lemma [32] for their mean values. Under the assumption that 
£,{X) e iJ^(R), the first of these vanishes up to order o(l) as E(/3) = 0, at least 
for P-a.e. realization w. What is clear is that the fluctuations of (|3.2|) will play an 
important role in the derivation of the appropriate Hamiltonian equations of motion. 
The second integral (|3.3p is less straightforward, as the mixing condition (|3.7p is 
in competition with the integral operators represented by the Fourier multiplier 
operators of the expression. We have that 

{P{x)D, i&nh{hD,)p{x)) 1^^^ \Dx^{X)\'' dX 

^ E(/3Z?,tanh(/iI?,)/3) j \Dx^{X)\'' dX. (3.16) 

There are two things to discuss with this statement. The first is that whenever 
7(a:,w) S is stationary with regard to some probability space (ri,7W, P), then 
an order zero Fourier multiplier operator applied to ^{x) is also stationary. Indeed, 
translation is respected 

m{D^)-i{x,TyLo) ^ ^ J e'''^''^'''^m{k)-f{x',TyLj)dx'dk (3.17) 

= — [ e'^^''-'-''^m{k)-i{x' -y,uj)dx'dk 
2tt J 

^ _L /" e"'(i^-y)-^')m{k)-/{x,u;)dx'dk 
27r J 

= m{D,)-f{x^y,Lo) . (3.18) 

Furthermore, continuous functions of 7 G C^, such as (7(7) — {'ym{Dx)j){0) are 
measurable. By the ergodic theorem, for any bounded measurable F 



lim J f F{T^g{^))dx^E{F{g)) 



and therefore the process Txg{"f) is ergodic. Secondly, the expectation values of 
quadratic functions of 7 may be computed from the covariance function p-y of the 
stationary process. For example, 

E(7to(L»j,)7) = lim E{j{x)m{D.j;)j{x - y)) 

= lim E{m{-Dy)j{x)-f{x - y)) = lim m{~-Dy)p^{y) 

= mi-Dy)py{0) . (3.19) 
Using these two facts, p.l6p is verified as the principal contribution from integral 

(ESI). 
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Lemma 3.3. [3] Suppose that (3(x;uj) is a stationary ergodic process which is mix- 
ing, with a rate a{y) which satisfies the condition (|3.7|) . Assume that E(/3) ~ and 
that (T/3 7^ 0. Define 

YM{X) = ^ [3{y)dy . (3.20) 



Jo 

As e tends to zero, we have, in the sense convergence in law that 

Y,{P)[X) ^ B{X) , (3.21) 
where B^{X) = B{X) is a normalized Brownian motion. 

In particular, let ,f{X) e 5 be a Schwartz class function, then 

Z,{p,f):= r -J—f3{-)f{X)dX (3.22) 



r'+oo /'+OC 

Yl{l3){X)f{X) dX - / -dxf{X)B{X) dX + o(l) 



This is to say that under the mild condition of mixing given in (|3.7p , the integrals in 
question converge to a canonical stationary process, for which only two parameters 
are distinguished, the mean value E(/3) and the variance cr^. This canonical process 
is given by white noise, 

f3{-)f{X)dX= / {£{p) + ^eapdxB{X))f{X)dX + o{^e), (3.23) 

£ J — oo 

where the equality is in the sense of convergence in law. The function f{X) in the 
integrand must be sufBciently smooth for the latter quantities to have a mathemat- 
ical sense. In fact we consider the operation of multiplication by /?(— ) to be in 
the distributional sense, which has for a limit the distribution y/eai3dxB{X) G S'. 
This is given a precise statement in the following lemma. 

Lemma 3.4. As a distribution, multiplication by j3{X/e) has a canonical limit in 
S' . Indeed, for f G S, 

P{^)f{X) = E(/3)/(X) + ^£a0dxB{X)f{X) + o{^e) . (3.24) 

Proof. Test the quantity above with a Schwartz class function ^p{X)] 

f3{-)f{XMX)dX (3.25) 

e 



= E(/3) J {fiXMX)) dX - ^eap J B{X)dx{fv) dX + o{^£) 
= f {E{p) + V^a^dxB{X))f{XMX)dX + o{V^) . 



This is to say that for each /, the random variable Z^{P,f) given in (|3.22|) is 
asymptotically normally distributed. Given two functions f,g€S, the covariance 
function E{Z^{P, f)Z^{P,g)) can be computed in the limit as e — > 0. Indeed 

E{Z,{(3,f)Z,{P,g)) = lJJpp{2£^)f{X)g{X')dXdX' 
= / / Ppix')f{X)g{X - ex') dXdx' 

^ J J Ppix')fiX){g{X) ~ ex'dxgiX) + -^x'^dlg{X) + ... )dXdx' . 
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Noting that the term at order e vanishes because pp is an even function, we have 
E(Z,(/3, /)Z,(/3, g)) = J pp{x') dx' J f{X)g{X) dX 

(3.26) 

J x'^pp{x') dx' J dxf{X)dxg{X) dX + ... 
In the hmit as e — > 0, this quantity converges to 

E(Zo(/)Zo(.g)) ^4 I f{X)g{X) dX , (3.27) 

where Zo{f) — J f{X)dxB{X) dX . This expression is consistent with the 

covariance of the white noise process being given by cr'^6{X — X'). □ 

In the case of a process I3{x) for which = 0, the hmit process for Yir{X) is of 
a different character. In particular, consider a stationary mixing process which is 
the derivative of another stationary process. Indeed let 7(x) € C"'+^(]R), and set 
/3(x) = 9^7(0;). Automatically E(/3) — and (jfj — 0. In this situation we have a 
different asymptotic result for the behavior of integrals such as in (|3.ip . 



Lemma 3.5. Suppose that 7(2;) G C""^"'^(R) is a stationary ergodic process which 
satisfies the mixing condition (|3.7p . and set P{x) = 9^7(2;). Then the process 
(3{X/e) is asymptotic in the sense of distributions to higher derivatives of Brownian 
motion. That is, for (p{X) € S we have 

PijMX) dX = e^'+i/V^ J d'x^'BiXMX) dX + o{e'+^^^) . (3.28) 

Proof. Using (p{X) as a test function, 

(3{-MX)dX - I d:^{-)^{X)dx 



e . £ 



= i-iye'- Jl{j)d^xV{X)dX 

= f Y,{-t){X)d^+\{X)dX 



= £'"+i/V^ d''+^B{X)^{X)dX + o{£'+^''') 



□ 

There are further technical results that we will use repeatedly in the analysis of 
the equations in the KdV asymptotic regime, having to do with limits in the sense 
of tempered distributions of products of scaled processes. In this context, consider 
7 = (71,72) a vector of stationary processes which satisfy the mixing conditions 
p.6p(|3.7p . Consider their product ')i{X/e)^2{[X+ct)/e) for some nonzero constant 
c as a tempered distribution in the limit e ^ . Define the covariance matrix of 
the vector process by 

C(7)=f^i 

\Pl2 O 

where 

/•oo 

a| = 2y^ E(7,(0)7,(y))d2/, Pi2 = P2i-/ E(7i(0)72(2/))dy . (3.29) 
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Lemma 3.6. If the vector process 7 = (71,72) is stationary and satisfies the mixing 
conditions (I3.6pp.7|) . then the process 

n(7) = 7i(y)rf2/,^° I2[v)dv) (3.30) 

converges to the two-dimensional Brownian motion B{X) = {Bi{X), B2{X)) with 
covariance matrix C(7). 

This result is analogous to Lemma 13.31 in the vector process case. From it, we 
derive the next useful result on products of two mixing processes. 

Lemma 3.7. Suppose that {f3i{x), P2{x)) is a C^{R) vector stationary ergodic pro- 
cess which satisfies the mixing condition (j3.6pp.7p . and let c be a nonzero constant. 
The new process formed by the product e~^Pi{X/s)P2{iX + ct)/£) converges in the 
sense of distributions on space-time to products of derivatives of a pair of Brow- 
nian motions with covariance matrix C(/3). More precisely, for a test function 
(p{X, t) e 5 then 

JPiif)P2{^Mx,t)dXdt 

= £/ dxBi{X)dxB2{X + ct)(p{X,t) dXdt + o{e), 

where the covariance matrix of {Bi{X), B2{X)) is given by C{f3). In case j3j = 
dx^j for indices j — 1,2, with jj g C^^'^^{M.) (so that a p. —Q if rj 7^ Oj the new 
process satisfies 

JPi{f)P2i^Mx,t)dXdt 

= £'■1+^2+1 J d'^+^Bi{X)d'^+^B2{X + ct)ip{X, t)dXdt + o(£'-i+'-2+i), 
where {Bi{X), B2{X)) are C{'y)- correlated. 

Proof. Start with the case in which both ap^ are nonzero, and write 

Jf3^{f)f32{^MX,t)dXdt ^Jf],{f)p2{^MX, ^) ^ 

= e'Jdx (/o^ /3i (r) dr) dx' (# /32(r') dr') ^{X, ^) ^ 



(3.32) 



■S{V~ej^ /3i(T)dr) P2{T')dT')dxdx^^{X, ^ 



X\ dXdX' 



The latter expression is a continuous function of the processes YeiP) = (Fe(/?i), yi(/32)) 
of equation (j3.20p , which itself converges in law to two-dimensional Brownian mo- 
tion with covariance matrix C(/3) as described by Donsker's invariance principle. 
Therefore the asymptotic expression for (|3.3ip is given by 

e j dxBi{X)dxB2{X + ct)ip{X,t)dXdt (3.33) 

where Bi{X) and B2{X) are two copies of Brownian motions with the correla- 
tion matrix C(/?). The general case reduces to the above particular case through 
integrations by parts. Indeed 

Pii-)P2{^^^Mx,t) dXdt 

e £ 

^r,+r. I dr^^^^^)dr^^^(^^±fL)^(^X,t)dXdt 
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(-1) 



Tri+r2 



71 ( f )72 ( ^{X, t) dXdt, 



which reduces the problem to the previous case. 



(3.34) 
□ 



There is another integral that needs to be evaluated in our further analysis. It 
has the form 



X + ^/ght ^ n 

P( — )I3{-Md, X, t)d9dXdt. 

e £ 



(3.35) 



I x,t J X 

The next lemma shows that such integrals have probability measures whose weak 
limits converge with order at least 0{e). 

Lemma 3.8. Suppose that {j3i{x) , (32{x)) is a C^{M.) vector stationary ergodic 
process which satisfy the mixing conditions (|3.6[) and (|3.7p . For test functions 

ipie,X,t)eS, 



dXdt 



X + -/ght 



X 



/3i(-)/32(-) + /32(-)/?i(-) 

£ e e e 



ip{0,X,t)d9 = 0{e). (3.36) 



Proof. The integral is written as the sum of two terms, each one of the form 

/ dxdt /^+^* A (f )/3, ilMo, X, t)de 

= £ J dXdt /^+^^* dx [V^jf A(s)ds) de (V^/o^ /3j (s)ds) ^{9, X, t)de 



:JdXdtJ^+^'' 



U I3^{s)ds) {V^ ij f3j{s)dsyxev{e,X,t)de 
V? A(s)dsj /o = l3j{s)dsjdgip{X + ^t,X,t) 
~(V^Jo^ PA^)ds)deviX,X,t) 



: / dXdt(^^ jf A(s)ds) [(V^ Jo^^^^ l3j{s)ds^dx^{X + ^t,X,t) 

-{V~e P,{s)dsyxv{X,X,t) 
JdXdt(y^jf (3,{s)ds) [^/3^{2£±^)^^x + V9ht,X,t) 

-^f3AfMx,x,t)_ 

(3.37) 

with i, j G {1, 2} and i ^ j. All of the terms have distributional limits which are 
at least 0(e). Simple cases which illustrate the estimate are: 

/: = eJdXdt[[^J^^ P^{s)ds)^P,{f) + [^J^'^ P2{.s)ds)^P,{f)]^{X,X, 
= I / dXdtdx[V^ lo^ I3i{s)ds^ jf l32{s)ds^v{X,X,t) 
= -f / dXdti^-^ if (3i{s)ds^ if l32{s)ds^dxv{X,X,i). 

(3.38) 
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//: = eJdXdt(^^jJ'p,{s)ds^^l3,i^^±^)ipiX + ^t,X,t) 

= -eJdXdti^^J^^(3,{s)dsj^dti^V^Io ' Pjis)dsjip{X + ^t,X,t) 

V^U Pds)dsji^V^Io ' P3is)dsjdMX + ^t,X,t). 

In these expressions, notice that the factors that appear are continuous functionals 
on path space. Therefore, as e i— > 0, they converge in law to functionals of Brownian 
motions pJj. The other remaining terms are easy to estimate. □ 

3.3. Random characteristic coordinates. Our method to derive the long-wave 
limit gives rise to a version of the KdV equation which has coefficients which are 
realization dependent. That is, the approximation process does not fully homoge- 
nize, and there are persistent, realization dependent effects that are as important 
as the classical effects of dispersion and of nonlinear interactions. The principal 
manifestation of this is the random overall wavespeed, expressed in the limit as 
£ ^ as 

3 /2 

co(X, Lo) = - t_^dxB(X) - e'^aKdv) ■ (3.39) 

The constant axdv is an adjustment to the characteristic velocity that is to be 
determined by an asymptotic analysis. The normally expected procedure is to 
solve the characteristic equations with this given wavespeed; 

^=co(X,co), X{Q) = Y , (3.40) 

to obtain characteristic coordinates (F, t) describing a net translational motion 
about which the more subtle nonlinear dispersive evolution takes place. In the 
context of a random bottom environment, however, the characteristic velocity field 
Co(X, cj) in p.39p has a component which is white noise, and when the flow of the 
characteristic vector field (|3.40p is required, (|3.39p is too singular to be able to make 
sense of a solution. 

Our derivation of the KdV equation is nonetheless performed in characteristic 
coordinates. To do this, our alternative strategy is to use a natural regularization of 
the characteristic wavespeed given in p.39p as an approximation, and to consider 
the characteristic coordinates indicated by ()3.40p to be the limit as e of a 
sequence of more regular flows. The regularized characteristic vector field that we 
use is 

dX , , r-r f £ ,X . , \ 
— = Ce{X,uj) := y/gh[l - — ) - e^aKdvj 

X{Q) = Y. 

We remark that as long as f3{x,uj) G C^(M) for P-a.e. realization u), the charac- 
teristic vector field Cg{X,uj) is C^, and for a given realization oj it is uniformly so 
in e. Therefore the flows X{t) = $f(F, w) exist for all e, and lie in a bounded 
subset of C^. The characteristics X{t) are themselves C^, and they are ordered by 
their initial values; if Yi < Y2 then for all t, Xi = cE>f(Yi,w) < X2 ^ ^'i{Y2,u}). 
As £ ^ there will normally not be a limit of the flows, but by standard 
compactness arguments there are limits 3'((F, w) in any C"(R), < a < 1 which 
converge uniformly on compact sets, and which preserve the ordered property of 
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the characteristics. Each such hmit X — <&t (F, uf) can be taken to be a weU-defined 
continuous and continuously invertible transformation. 

To understand the asymptotic behavior of the transformation to characteristic 
coordinates, write 

^{Y, Lo) = + eX^{t) + e''X^{t) + . . . , (3.41) 

where X^{Q) = Y , and X\{)) = for j > 1 provide the initial conditions for the 
flow. Substituting this into the characteristic equation gives the result that 

= V^, X\t) = Y + t^, (3.42) 

at 

dX\t) _ 1 /7^,^o 



dt 2 \ h 



-Ji^P(X"{t)/s,co) 



= '\J{m + t^/9h)le,Lo), (3.43) 



thus, variations to the characteristics are given by 

,(Y+t^)/e 



^'W = -^/ Ks,u;)ds. (3.44) 



lY/e 

The final term relevant to our considerations is 

dX^{t) 



= -ygh O'Kdv (3.45) 



dt 

which integrates simply to X'^{t) — —\fghaKdv t. Studying the integral expressions 



for X^{t) more closely, we find that 



^'(t) = -^( / mds- I mds), (3.^ 



which converges in law to Brownian motion as e ^ 0, according to our discussion 
in Section [321 Hence 

X\t) = -^{B{Y + ty^)-B{Y))+o{V^) . (3.47) 

In particular, the term eX^(t) contributes at order e^/^. Due to Brownian scaling 
and to the property of independence of increments, 

X\t) = -^B^^y^it^) = -Vi(f^^)i?c.(y)(i) . (3.48) 

We note that the realizations lo{Y) of Brownian motion depend on the different 
initial positions Y, and in particular that for distinct initial points Yi and Y2 the 
selection of realizations -Bu;(i'i)(0 ^^'^ ^w(y2)(^) Brownian motion are indepen- 
dent, as long as Y2 — Yi > y/ght. Putting this information together, an expression 
for the characteristic flow is given by 



^■z AY+tVgh)/e 

X{t,Y-e,Lo) = Y + t^-— P{s,Lo)ds 

^n- JYje 

-^fghaKdve^i + -- - ■ (3.49) 
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As £ tends to 0, the characteristics tend to the hmiting distribution of paths given 

by 

3 /2 

X{t,Y-uj)^Y + t^~^-^^^B^^Y){t)-e^V9haKdvt + --- ■ (3.50) 
Inverting the expression gives a formula for Y in terms of X and <; 

Y{t,X-e,Lo) = '^U{X;uj)^X-t^+— / P{s,uj)ds 

+e^^/ghaKdvt-\ . (3.51) 

As £ tends to 0, 



3 /2 

Y{t,X;Lu)^X -t^+^^^^B^(^x){t) + V9haKdveh + -- - . 
The Jacobian of the flow has the fohowing asymptotic expansion 

In the Hmit as e tends to zero, the Jacobian (|3.52p . when muhiplying a test function, 
behaves asymptoticaUy as 



dX , e^/^aa 



^-'-^VghdxB^^yJt). (3.53) 



4. BOUSSINESQ REGIME 

We now return to the expression (|2.23p for the scaled Hamiltonian, in order to 
give a formal derivation of the appropriate Boussinesq system in this regime. Recall 
that /3 is a mean zero, stationary mixing process with correlation function pfj. Using 
the analysis of the previous section, we write the leading order contributions of the 
second and fourth terms of (|2.23|) in the form 

Pi-,u;)\DxaX)\^dX =^£ap / dxB{X)\Dxi{X)\^ + o{^e) (4.1) 



and 

{[3DA^iMhD,)[3){-)\Dxax)\^ dX = Eif3D^tanh{hD,)P) / \Dx^{X)\^ dX , 

-OO ^ J 

(4.2) 

where as in p.l9[) we calculate that 

£{(3D^ta.nh{hD^)l3) = {Dyi&nh{hDy)pp){0) ap. (4.3) 

The constant will contribute to an adjustment of the linear wavespeed in the 
Boussinesq regime. The equalities in (|4.1I) (|4.2|) are to be taken in the sense of law 
of the corresponding random processes. 

We now include these expansions into the Hamiltonian, retaining terms up to 
order 0{e^) and dropping those of higher order 

' ({h - e'"^npdxB{X) - e''ap)\DxS,\'' + gv^)dX 



J {iDxvDxi - Y^DU)dX + o(£^). (4.4) 
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We note the role of a stochastic efFective depth played by 

ho{X) e'^''^(TiidxB{X) - e^ap + o{e^) (4.5) 

which is a function of the long length scale variables alone. Since it is normally not 
necessary to introduce characteristic coordinates in this derivation, a regularization 
such as described in section 13.31 is not required, and the limiting effective depth 
/io(X) is used directly in the averaged Hamiltonian. 

Changing the variables {r],£) to (77,1* = dxO^ the Hamiltonian becomes 

Hi = '-^J {hoiX)u'+g7^^^e\^idxuf-w^))dX. (4.6) 

The symplectic structure has to be modified accordingly as in [5]. Consider the 
transformation w v = f(w), which transforms Hamilton's equations 



to the form 



dtw = J5^H{w) (4.7) 
dtv = Ji5,Hi{v) (4.8) 



with a new symplectic structure 

Ji = d^fJid^ff , (4.9) 



where d^f is the Jacobian of the map /. In our case, w 



„ I uy, and the matrix J = e , „ is transformed to Ji = £ 1 „ 
dx J \~I OJ \~dx 

where the power of e is due to the scaling transformations (|2.15p . The evolution 

equations take the form 



In the end we find the Boussinesq system in the form 



3 

dtu — —gdx'H — £^udxu . (4-11) 



dtV - -dx{{K{X)+e''ri)u)-e'^—dlu, 



While this form of Boussinesq system appears most naturally from a direct ex- 
pansion of the Hamiltonian of the problem of water waves, the resulting system 
of partial differential equations is not well posed, and it is rarely used directly in 
modeling. In the present setting, the situation is further aggravated by the fact that 
a coefficient in the above system is singular, as it involves the second derivative of 
a Brownian motion. Several routes to resolving these issues are possible, modifying 
the linear dispersion relation for the Boussinesq system, and regularizing the coef- 
ficients as in section 13. 3[ for example. However we will not pursue this direction 
of inquiry in the present paper, preferring to make a more systematic study of the 
KdV scaling regime. 
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5. The KdV regime 



In the case of the Boussinesq derivation, the hmit of certain integrals in the 
water waves Hamiltonian will give rise to singular coefficients in the resulting equa- 
tions of motion. This is even more true in the case of the KdV regime; indeed the 
transformation to characteristic coordinates will give rise to a modified symplectic 
structure which involves a second derivative of Brownian motion, something that is 
not acceptable on an analytic level. To get around this difficulty, we regularize the 
linear wavespeed as described in section [3.31 a process which consists of retaining 
certain terms with rapidly varying coefficients in the Hamiltonian, and only taking 
the limit after the long wave equations are derived. We assume that ap > 0, which 
implies that the resulting realization dependent fluctuations are maximally signifi- 
cant in the limit, and we will perform the smoothing procedure in a way which is 
consistent with this assumption. 



5.1. Successive changes of variables. We start again from the expression (|2.23p 
for the Hamiltonian. As in the derivation of the Boussinesq system, we first change 
the variables (ry,^) to (ry,M = dxS,), leading to a transformed Hamiltonian Hf 
defined by 



-dx 



dX 



and a modified symplectic structure Ji = e 
variables is defined by the transformation 



-(r + s) 



-dx 



-(r-s) 



(5.1) 

The next change of 
(5.2) 



The new symplectic structure resulting from this transformation is 

1 dxh, 

I -Ox 

-3 



J2 



4 



4 



d 



(5.3) 



'X 



whose off-diagonal terms quantify the scattering of solutions due to variations in 
the topography. In this expression, we retain the regularized expression : 



e/3( — ) - 



(5.4) 



The Hamiltonian is written as 

-3 



r^s - rs'^ + s^)^dX + o{e^) 




(5.5) 



Notice that, except for the first term in the Hamiltonian, he appears in terms that 
are already of order and thus can be replaced there by the constant h in this 
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asymptotic calculation. Denoting by 
we rewrite i/f in the form: 

-3 



C2 



2 V 4/i' 



^eir^ + s^) - c,s^[{dxrf - 2{dxr)idxs) + (Oxs) 
+026^ (r^ - r^s - rs^ + s^^ dX + o(e^) . 
Hamilton's equations for (r, s) take the form 



dt 



J2 



5sm 



where J^-fff ^^-fff ^-^^ computed as follows: 

5^H^ = 4(VFv2r + cie2(2a|r-2o»|,s)+C2e2(3r2-2rs-s2) 

5si?2 = 4(\/5^2s-ci£2(2a|r-29|s)-C2e^(r2+2rs-3s2) 
Hamilton's equations are explicitly 



(5.6) 



(5.7) 



(5.8) 



dtT — 



X 



2rs 



1 dxhe 

4 h. 



+ {ci{d\s ~ d\r) + \c2{~r'^ ~ 2rs + 3s^] 



dfS = 



fghls + e'^{ci{d\s - d\r) + \c2(~r'^ - 2rs + Ss^) 



(5.9) 



1 dxh, 
4 /i. 



fghlr + £2 (^ci(a|r - d\s) + ic2(3r2 - 2rs - s^)^ 

(5.10) 

In the action of J2SH2, there are products of and its derivatives, and each 
factor tends to a distribution (see Lemma [3. 2p in the limit e ^ 0. The product is 
nevertheless well defined because of the form it takes: 

J dxK h,-^'^f{X)dX = 2 j dxhl/^f{X)dX . 



We perform an additional change of scale of s relative to r defined by 



1 

£-3/2 



(5.11) 



which puts forward r{X,t) as the main component of the solution which is an- 
ticipated to be traveling principally to the right, with a relatively small scattered 
component si{X^t) propagating principally to the left. The transformation leads 
to a modified symplectic structure 



1 



-dx 



1 dxhe 



4£3/2 

1 dxhe 1 „ 

-^ox 



(5.12) 



4£3/2 he 
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and a final Hamiltonian 



(5.13) 

The equations stemming from the Hamiltonian (|5.13p and the above symplectic 
structure are 

dtr — ~dx s/ghlr + {cid\r + ^C2p) 



+e^(— £2Ci9^si — £2C2rsi — \e'^C2s\] 



'4 



(5.14) 



+C2( 

dx 



2' 



.V\s^ + 3^3^2)) 



+ C2( 2 

1 ax/3(f ) 



■rir'^ — rsi + ^eis\) 
e^^y/gh^r + (ci(9|-r — e^dxSi) 
+C2(|r^ - eirsi - y^sD) 



(5.15) 



It is ambiguous at this point precisely which terms of the above system of partial 
differential equations play a role in the asymptotic description of solutions in the 
limit as e tends to zero. The transformation (|5.1ip is not homogeneous in the 
perturbation parameter £, and because of fluctuations there are numerous cancel- 
lations that occur in the remaining terms, not all of them having an influence on 
the asymptotic regime (see Lemmas 13.21 and 13.51 for example) . We will show in 
the subsequent analysis of Section 15.31 that the asymptotic behavior of solutions 
of equations ()5.14p ()5.15p as £ —> is governed by the following coupled system of 
equations, with an appropriate choice of the parameters axdv and b. 



dtr 



-dx 



e br 



(5.16) 



dtSi = ^/ghi 



ghdxs, + -.[^e-^'^d^(3{-)r, 
4 V ft £ 



(5.17) 



where the regularized velocity is Ce(X) = y/gh{l — ■^P{X/e) — e'^axdv)- There are 
two free parameters in this system of equations, namely, aKdv and b. They will be 
determined by the consistency analysis of Section [5.31 as fixed points of the solution 
process and the asymptotic analysis. In the end we find that 



CKdV 



1 

2h 
7ci 
'64/i3 



E((9./3)^). 



(5.18) 
(5.19) 
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5.2. Solution procedure for the random KdV equations. In this section 
we describe a reduction procedure for the system of equations (|5.16p - ()5.17|l that 
expresses the solution component r(X,t) in terms of a solution q{Y,T) of a deter- 
ministic equation similar to the KdV equation, under a random change of variables 
{Y X{t,Y)) and a scaling r = e'^t to the KdV time. The scattered component 
si{X,t) is an expression involving integrations along characteristics. The solution 
depends upon the two parameters andv and b. We retain the regularized form of 
the characteristic velocity c^{X), only taking the limit as e ^ in expressions for 
the solution. 

Substitute r — dxR into (|5.16|) : the resulting equation for R is 

dtR = -Ce{X)dxR " e^icid^R + ^c2{dxR)^) + e^bR. (5.20) 
Transform to characteristic coordinates as in Section [3.31 

^ = c,(X), XiO) = Y. (5.21) 

We denote the flow hy X = <i>f(y), which is a regularized realization dependent 
change of variables. Define Q{Y,t) = R{X,t) so that Q satisfies 

drQ = -cid^Q - ^C2{dYQf + bQ . (5.22) 

To solve the initial value problem, set q{Y,Q) = r{Y,Q) = r'-\Y), and solve the 
deterministic equation 

drQ = —cidyQ — 3c2qdYq + bq (5.23) 

for q{Y,T) = dyQiY^r). If 6 = 0, equation (|5.23p is the classical KdV equation. 
Additionally, for each realization /3(x,w) the regularized ODE (|5.2ip defining the 
flow has a solution given hy X — X{t,Y]£,Lo). With these two ingredients, the 
solution r{X,t) of equation (|5.16p is given by 

r(X, t) = dxQ{Y{t, X- e, uj), eh) = dYQ{Y{t, X; e, w), eh)dxY{t, X; e, u;) (5.24) 

where 9yX(i, F; e, w) is the Jacobian of the flow (|5.2ip as described in section 
13. 3i and dxY(t, X;e,LLj) is its inverse. This is an expression of the solution of the 
regularized equation. 

The equation (|5.17p describes the scattered component of the KdV system above, 
whose solution is expressed by integration of a forcing term which is given in terms 
of r{X, t) along left-moving characteristics . Explicitly, 

si(X,t) = s\{X + ^t) 

+^yi/o 9,/?( ^+^(*-'') )r(X + ^{t - i'), t') dt' (5.25) 
= .;(X + V^i) + 4/^+^*5./3(f)r(e,t+^)d^?. 

The small parameter e is still present in the regularization; to complete the descrip- 
tion we consider the limit of the expressions (|5.24p(|5.25p as e tends to zero. The 
solution of (|5.22p is smooth, and admits a Taylor expansion in its arguments. The 
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inverse Jacobian has an asymptotic expression as well. Therefore, one writes 
riX,t) = dxQ{Y{X,t-uj),eH) = dYQ{Y{X,t-uj),t)dxY(X,t-oj) 

= q(X V¥^t,eH){l + 4(/3(f ) - I3{^^) 

= q{X -^t,e^t) 

+dx (q{X - V^t, sH)[f, /(L^,)/, m dt') )+0{s'). 

(5.26) 

Proposition 5.1. In the limit as e tends to zero, the expression (j5.26p for the 

solution of (j5.16p is asymptotic as a distribution to 

r{X,t) = q{X ~ yfipit.e^t) 

3/2 / \ , (5.27) 

V-^<f^dx[q{X - V9ht,eH)B^^x){t)) + o{e^/^). 
The expression for (j5.25p for the solution si is asymptotic as a distribution to 
si{X,t) = s\{X + ^t) 



Aha 



X 



B{e)-^q{20-X 



X — 6 
ght,e'^{t-\ ^=^))d9 



gh 



4h. 



^ (^dxB{X + V9ht)q{X + Vght, 0) - dxB{X)q{X - Vght, eH) 

-2K^{b{X + ^t)dxq{X + V^t, 0) - B{X)dxq{X - ^t, sH) 

(5.28) 

Proof. The expression for the limit of ri follows directly from the application of 
Lemma 13.51 It is an expression which exhibits both randomness in its amplitude, 
as well as in location as per the random characteristic coordinates in which it is 
expressed. For the calculation for the limit of si, we substitute the expression (|5.26p 
in ■■ 

si{x,t)^s'i{x + Vght) 

x+^/^t 



^/ d.,f3C-)q{2e-X -^t,e^t 



x-e 



X 

X + ^/ght 



))de 



X 



I3{s)ds 



(5.29) 



'.dxq{2e- X -^/^t,e^{t 



de 



_l X + y/ght 
£ 2 " 



X 



' 2 ^s' 



29 -X 



ght 



xq(29 - X - y/^t,e^{t 



X- 



))de. 



Except for the first term s\ that remains unchanged, all the terms appearing in 
the limiting expression (j5.28p come from the first integral in the expression of si, 
where we performed several integrations by parts and use the fact that dtq{X,e^t) 
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is 0{e'^). By more integration by parts, using the fact that dx = sdx we can show 
that the third term (third and fourth hnes) in the expression (|5.29p is ©(e^/^). 
Let us turn to the last term (fifth and sixth hnes) of (|5.29p . For the term con- 
taining dx0^{O / s), integration by parts wiU produce an additional e and the term 
will eventually be of order ©(e^/^). To estimate the term containing the product 

dxP{j)P{ ^^~^~"^* ), the integration by parts moves the derivative dx to all other 
terms. The only contribution that will not produce an e is when the derivative acts 



d^Pi ) = i )• (5-30) 



After some simple manipulations, this term is again 0(e^/^). □ 



5.3. Consistency of the resulting system of equations. In this subsection, we 
complete the cycle of a self-consistency analysis for equations (|5.16p and (|5.17p . out 
of which the two so-far undetermined constants axdv and b are selected. It is clear 
that not all terms in equations (j5.14p and (|5.15p are of equal importance in the limit 
as e ^ 0. Recall the criterion as presented in section [3l which states that a term 
a{X,t;e,uj) is of order 0{e^) if for any space-time test function (f{X,t) £ S the 
measures Pe induced by J a{X,t;e,u!)(p{X,t) dXdt converge weakly to a limit 
Pq as e tends to zero. In the present case, the analysis consists of (i) the derivation 
of an expression for the solutions of (|5.16p(|5.17p which are stated in (|5.25p and 
(|5.26p . and depend upon the two parameters axdv smd 6; (ii) the examination of 
the terms in (|5.14p . including in particular those which do not appear in (|5.16p (re- 
spectively, all the terms in (|5.15p . in particular those that do not appear in (|5.17p ). 
Using the expressions (|5.25p(|5.26p we then show that, except terms which appear 
in (|5.16p (respectively (|5.17p ). they are asymptotically of order o(e^) (respectively, 
of order o(l)). Both the system (|5.16p(|5.17p and the solution expressions ()5.25p 
()5.26p depend upon parameters aKdv and b. (iii) The demonstration that these 
constants can be chosen so that there is a fixed point of this analysis. Namely, the 
solution depending upon the constants a^dv and b has asymptotic behavior which 
satisfies the equations (j5.16p - (j5.17p with the same choice of constants. 

Let us denote the terms in ()5.14l) by 

Ir = e'^dx(^ - e'-^/'^cidx-si - e^/^C2rsi - ie^C2Si^ 

III, = _i£^e2ci(-air + e3/2a|,^) 
IVr = -I^4(±)e2c2(-ir2-e3/2rsi + |e3s2). 



(5.31) 



Similarly, we denote the terms in (|5.15p by 
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Is 



II, = i^£-3/2^, 

III, = 1 ^^£1/2^,(92^^ _£3/2g|,^) 

IV, = ^^ille^^c^dr^-s^/^rs.-^e^sl). 

(5.32) 

The purpose is to evaluate the asymptotic behavior of each of these terms as e ^ 0. 
Lemma 5.2. The term II,- has the asymptotic behavior 



llr = —^le'£{P')dxr{X,t)+o{e'). (5.33) 
Lemma 5.3. The term 11, has the behavior 

^^s^'-^^|l^.^3{^)r + o{l), (5.34) 
and this expression has an asymptotic limit as e Q which is 

\^af,dlB{X,u)q{X-^t,T). (5.35) 

Lemma 5.4. 

III. = ^e^E{{d.pf)dxr - ^e^£{{d.(3f)r + o{s'), (5.36) 

III, ^ e-^^HUr ^ 0(6^/2). (5.37) 
Lemma 5.5. The remaining terms have the following asymptotic behavior 

= o{e^) , IVr = o{e^) (5.38) 

and 

Is = o(l) , IV, = o(l). (5.39) 

Lemma 5.6. Finally the linear term —dx{^/ gh^r) in the equation f5.16\) has the 
asymptotic behavior 

- dxiVghsr) = -^dx [(l - |^/3(|) - + ^E(/32)))r] + 0(6^). (5.40) 

The proofs of these lemmas are the content of Section [5^ Using these asymptotic 
results in system (|5.14l)(|5.15l) . and retaining only the leading terms, it reduces to 
(|5.16p(|5.17p . with possibly different parameter values. When the parameters are 
chosen appropriately, the asymptotic behavior of the equations matches that of the 
solutions and the consistency procedure is closed. 

Theorem 5.7. The result of the consistency analysis is that the free parameters in 
equations i5.16])(KJ7\ ) are 

aw ^ ^a. + ^E(/5^) + ^E((a./^)^), (5.41) 
E((9,/?)3). (5.42) 



64/i3 
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The parameter aKdv represents an adjustment at ©(e^) to the overall wavespeed, 
while the sign of b governs the stability of solutions. In many cases, h vanishes. 

Proposition 5.8. If the statistics of the ensemble (ri,A^,P) are reversible in x, 
then 6 = 0. 

By reversible, we mean that the inversion x —x preserves the probability 
measure P, implying that £((9^/3)^) = 0. 

5.4. Proofs of the above lemmas. In the analysis of the numerous integrals 
that go in to this consistency result, it is convenient to use the bracket notation as 
shorthand for integrations; 

(/,5> := / f JiX,t)g{X,t)dXdt . 

Proof of Lemma \5.Si We first rewrite 11^ as 

II. = ^e^f^dx{Vh'e)s, = ^ey^dx{Vh'e-E{Vh'e))si. (5.43) 

For any test function ip{X,t), we compute (ip^llr) by substituting the expression 
(|5.25p for si. This gives two terms, the first being 

£3/2^(^,9^(7^- E(V^))s;) =-^e3/2^(y^_E(v/X:)),ax(s?¥>)). (5.44) 
Since 

and because 

y/K- E(v/^) = Vh{l - ^Kj)) - Vh + 0{e^) = --£=/3(|) +0(£2), (5.45) 
the first term in ((^,11.) is of order o{s'^). The second term in the expression of 

(((9, Ilr) is 

t^(^,5x(v^-E(v/^)) / d,p{-)r{0,t+^^)de}. (5.46) 

By integration by parts, 

^ = (dx^^ iVK - EiVK)) /^+^* d,(3{l)r{9, t + ^)de) 

-||(^, {VK - EiVK)) ^ d.^p{^d,r{e^ t + ^)de) 

{Vh; - E{vr,)) [d,p{^£±^y{x + v^t) - d.,(3{§)r{x, t)] ) 
= ii9x - ^dt)vAv^e-EiVK))Ix^'^'d.0{l)r{e,t + ^)de) 

+ || (vp,(V^-E(V7i;))9./3(f)r(X,<)) 

(5.47) 

Analyze the second term first, 

1 /S",/,. n /;32n 



WATER WAVES OVER RANDOM TOPOGRAPHY 



27 



Replacing r by its expression ()5.26p 



^yie {^,dAP') e^dxqn 



+eq (/?(f)-/3(^^)) )+0(£f) 



(5.49) 



By integration by parts, the first term of A2 is o(e^). The second to the last term 
of A2 can be rewritten as 



(5.50) 



which again by integration by parts contributes to o(e^). The last term of A2 
contributes only o(£^) due to Lemma [3.71 Now turn to Ai. 

' X \ c-X + ^/ght oi e\ d 



^1 = -^.vi^Hidx ^,d,)^,miV''mfAo,t+^)do) 

+ lkVi^''ii9x - ^d,)^, /3(f) [/3(^±^)rO(X + ^t) - /3(f )r(X, t<\ ) 
+o(£2). 

(5.51) 

The first term in the second line of A\ is o(e^) due to Lemma [3.71 The last term 
of A\ is 

-1^ V^E(/32) ((5^ - ^a.)^, <? ) + o{e^) 
leading to a contribution to Ilr of 



-^|E(/J^)a..(^,0+o(e^). 
We now turn to the first term of A\ which we denote ^3, and write it as 

A, = ~^./f ((ax - vk^*)^'/5(f )/_^+^*/3(f)(ax 



We express (dx ~ '^^'^t)'^ terms of q as 



(5.53) 



(5.54) 



(dx - :^dt)riX, t) = 2dxq + ^ (a./3(|) - 2d,Pi ^ )) q + 0(e). (5.55) 

Substitution of the above in ^3 gives rise to three terms, (z), (ii), and (iii) which 
have the form (after we have dropped the constants): 

= ^'ii^x - -^dt)v, /3(f) /^+^* P09xq{2e -x-^t, sHt + ^))de) 



The term («) is of the form 



X + ^fght n 

Pi-)(3i-)^ie,x,t)de). 



(5.56) 
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Applying Lcmma lSTSt we show that this term is 0{e^) , and thus docs not contribute 
to the hmit of 11^. By integration by parts, the term (m) is 0{e^). FinaUy, for term 
{Hi), we write d^p C"-^;^* ) = --^ ^^( ^e-x-v^t .^^ leading to (Hi) being 
again of order 0{e^). 



□ 



Proof of Lemma[EM- Using that h^^h- efj{f) + 0{e'^) 

Since r{X, t) = q{X - ^/ght, eH) + 0(e), the second term of (|5.57p is 

'-^Jld,f3\^){q{X-^t,eh) + 0{e)) =0{e'/') (5.58) 



due to Lemma 13.51 Compute the limit as e ^ of 11^ . Substituting the expression 
()5.26p for r, we get, for any test function ip{x, t) 



3 

11. IT 

4 V h 



(^,9./3(f)^g(/3(f)-/3(^^))). 

(5.59) 

The first term of the RHS of (|5.59p tends to the first term of (|5.35l) by application 
of Lemma [3.51 The second term of (j5.59p is rewritten, by integration by parts, as 



Clearly all terms are o(l). The third term of (|5.59p is rewritten 

^Vf (^,9./32(f )9) - #71 (^,/3(f )5.^(M*),), (5.60) 
which is o(l) by application of Lemmas 13.51 and 13.71 



□ 



Proof of Lemma \5.4\ Decompose 111^ as the sum of the two terms 

C = — £ — J^^^xr, (5.61) 

We compute dj^r from (|5.26p and do not write terms that will clearly give a con- 
tribution of o(e^). We get 

C) = ^e^^, M(3a^qa./3 + e-\dl(3)) + o{e^). (5.63) 



The first term in (|5.63p . denoted Ci is 



Ci) = ^e^E{{d,(3f){^, dxq{X - ^t, eh)). (5.64) 
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We substitute dxq in terms of r in Ci using ()5.26|) and we write 



0{s). 



We then conclude that Ci can be written as a functional of r as 
The second term in C, denoted by C2 is 



(5.65) 



(5.66) 



L6/i 

Cl 



^,a,(a./3)2(l + -/3) g)+0(e3) 



16/i3 

We conclude that the term C of Illr is 

C=^^e^E{{d.,Pf)dxr- 



e^£{{d.,[3f){^,q)+o{e^). 



^e^E((a./3)3)r + o(s^). 



(5.67) 



(5.68) 



We compute the term D of III^ given in (j5.62p . For this, we compute d\si in terms 
of r and get: 

dls,{X, t) = dls\{X + ^t) + ^Le-I [ia2^(2£±YSM)^0(^ + 

+9./3(^M*)9^^0(^ + VF^t) - ldlP{^)r{X,t) - d^p{f)dxr{X,t) 
+ ^d^f3{^^±^)dtr{X + V^t,0)-^d^f3{f)dtr{X,t) 



+ jElx^'^''mi)9ur{0,t+^)d9 

(5.69) 

All terms containing the process /3 or its derivatives at two different points X/e 
and {X + y/ght)/e will not contribute because of Lemma l3.7l The term containing 
s° will be o(e^). The remaining terms that need attention are 



^(^, ^ [s-'dlP r{X, t) + d.,f3{dxr[X, t) + -^dtr[X, t))\ ) 
d^,'-eix''^''mi)dur{9,t+^)de). 



(5.70) 



cie" 1 
16 'gfP 



Noting that dxr + -^^dtr = 0{e), we have that 



cie 

'mi 



^ (^d^l3{dxr{X, t) + -^dtr{X, t)))) ^ 0{e'). (5.71) 



The first term in (|5.70p has the form 



32/1 



32P 



.(^, 5.((a,/3)2 - E((a,/3)2))r) - E((a,/3)3)r> + 0(^2) 



(5.72) 
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Integrating by parts the first term of ()5.72|) . we get two contributions; wlicn tlie 
derivative acts on (p, it is o(e^) using Lemma [3.51 and the fact that r = q + 0{e). 
When the derivative acts on r, we get: 

- ^{ ^, {{d.pr - Empr))dxr ). (5.73) 

Here we replace dxr by its expression in terms of q: 

dxr = dxq + ^l{9^f3i§) - d,P{ ^~^* )) + 0{e). (5.74) 
The resulting contribution for (|5.73p is 

- 'P, {{d.Pr - E{{d.f3r))dxr ) = (p, E((9,/3)3)r ) + o{e'). (5.75) 

The last term to consider is the fourth term of (I5.70p where the derivatives with 
respect to t can be moved outside the integral using the fact that 



^d^f3{^^^^)dtr{X + VPt, 0) - ghd^p{^^±^)dxr{X + ^t, 0) 
-f 9../3(^±^)r(X + Vght, 0). 

(5.76) 

Using Lemma 13.71 again . 



(5.77) 

Using the derivative in the first factor of dx(3 appearing in the above expression 
and integrating by parts leads to the appearance of an additional s, making the 
expression 0{e'^). We have obtained that 

D^~-^Emi3r)r + o{e'). (5.78) 

Adding the expression for C and D , we have shown that (|5.36p describes the 
asymptotic behavior of Illr. □ 



Proof of Lemma \5.5[ Following the criterion of Section 3, these terms are integrated 
against test functions f , and derivatives can be moved to ip by integration by parts. 

□ 

Proof of Lemma [57S[ The regularized depth is defined as h^{X) = h~ £/3(^) — 
e'^ap. Thus the regularized linear wave speed is 

V9^=v9;(l-|;^-.^|-e^^)+*"). (5.7!,) 

The term {ip.dxiP'^r)) is calculated as 

dxiP^r)) = £{p^)dxr) - {dx^, {p^ - E(/32))r). (5.80) 
Using that r — q + 0{e), we get that the second term in (|5.80p is 0{^/£). □ 
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6. Remarks on the expectation of solutions 

It is normal to calculate E{r{X, t,u!)) — p{X, t) as a basic prediction of the solu- 
tion r{X,t,u!) itself. We remark that r{X,t,Lu) is a realization dependent function 
where the randomness manifests itself on the same level as dispersive and nonlinear 
effects. In the paper 21 on apparent diffusion, the authors present an analysis 
of the function p{X, t) in the case of the linear water wave problem with bottom 
given by {y = —h + ^f3{X/ e)}. In the fully nonlinear regime of the present paper, 
diffusion is weaker, and occurs only on time scales larger than those of 0(1) in KdV 
time r, as the following calculation shows. 

In the sense of weak limits of probability measures, as e — > 0, 

r{X,t)=q{Y,T), (6.1) 

where 

g3/2 

Y = X - y^t+ —{ghy/^<Tf)B{t) +e^aKdv\/ght, and r = e^i. (6.2) 
Compute the expectation of the main component of the solution r : 

l-oo ^3/2 

E{r{X,t)) = I q{X - y/ght + ^^af3{ghY/'^u + e'^aKdv\/^t,T)diiB(t){'U') 

g3/2 ^ ^2 

q{X - \/ght + -^ap{gh)^u + e^axdv^/ght, T)e~'^du. 




(6.3) 



Assuming that max^ |(7(.,t)|/^i < c», we have for fixed t, 

1 f"^ £3/2 ^ 

maxE(r(X,i)) <max— == / \q{X' + —-afi{gh)^u,T)\du 



<'^hLl{gh)-i I \q{v,T)\dv. (6.4) 



This time decay of order e~^^^t~^/^ = (£r)^^/^ shows that the diffusion coefficient 
is of order 0{e), meaning that diffusion effects occur at an order higher that the 
one considered for the derivation of the KdV equation. To observe diffusion created 
by random effect at the order of the relevant terms for the KdV would require a 
scaling for the bottom variations of the form —h + ^/ef3{x,uj), which is a 'rougher' 
bottom that the one considered in this paper. This is the natural scaling that was 
considered in the linear analysis of [21]. However, such a hypothesis also affects the 
nonlinear and dispersive nature of solutions and indeed it will introduce additional 
terms in the nonlinear coupled system of equations for (r, s) that would have to be 
taken into account. This is beyond the scope of the present paper and is planned 
as the focus of a subsequent study. 
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